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Abstract 


Computational  fluid  d>aiamics  calculations  have  been  performed  for  a  multibody  system 
consisting  of  a  main  projectile  and  three  sabot  components.  Numerical  flow  field  computations 
have  been  made  for  various  orientations  and  locations  of  sabots  using  an  unsteady,  zonal  Navier- 
Stokes  code  and  the  Chimera  composite  grid  discretization  technique  at  H.  -  4.0  and  a  =  0°. 
Computational  grids  have  been  obtained  for  the  projectile  and  sabot  independently  and  then 
overset  to  form  the  complete  grid  system.  Comput^  results  have  been  obtained  for  sabot  angles 
of  attack  of  5, 10, 15,  and  25°.  Computed  results  show  the  details  of  the  expected  flow  field 
features  including  the  shock  interactions.  Both  laminar  and  turbulent  computations  for  the  25° 
case  predict  similar  results.  Computed  results  for  oth^  sabot  positions  are  compared  with  the 
experimental  data  obtained  in  Canada  for  the  same  configuration  and  conditions  and  are 
generally  found  to  be  in  good  agreement  with  the  data. 
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1.  Introduction 


During  the  gun-launch  process  of  modem  fin-stabilized  kinetic  energy  projectiles,  sabots  are 
utilized  to  reduce  in-bore  balloting  and  smoothly  carry  the  projectiles  in  the  gun  tube.  Once  free  of 
the  gun  tube,  the  sabot  components  must  be  discarded  to  reduce  the  drag  of  the  round.  Good  sabot 
separation  is  important  to  obtain  a  repeatable  launch  and  flight  of  the  projectiles.  It  has  been 
demonstrated  [1]  that  aerodynamic  interference  during  the  launch  process  can  adversely  affect  the 
projectile  trajectory  and  increase  on-target  dispersion.  Mechanical  interference  during  the  sabot 
separation  can  also  alter  the  projectile  trajectory  and  may  lead  to  unacceptable  loss  of  accuracy  at 
the  target.  The  aerodynamic  interference  of  the  projectile  and  sabot  flow  field  is  quite  complex  (see 
Figure  1)  and  involves  three-dimensional  (3-D)  shock/boundary  layer  interactions  and  separated  flow 
regions.  The  sabot  separation  in  actual  flights  can  be  asymmetrical  which  can  further  magnify  the 
interference  effects. 


Figure  1.  Spark  Shadowgraph  for  a  Projectile-Sabot  System. 
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The  sabot  discard  aerodynamics  have  been  studied  both  theoretically  and  experimentally  over 
the  past  few  decades.  Initial  analytical  modeling  efforts  [2]  were  empirical  and  based  on  Newtonian 
flow  approximations.  Although  not  general,  this  technique  accurately  modeled  the  flow  in  some 
cases.  A  more  recent  extension  [3]  of  this  modeling  effort  included  an  integrated  flow  element 
approach  that  utilized  local  shock/expansion  procedures  based  on  actual  experimentally  measured 
sabot  surface  pressures.  Again,  these  analytical  modeling  techniques  have  their  limitations,  and  in 
general,  these  theories  cannot  completely  represent  the  complex  flow  fields  associated  with  the 
interaction  of  the  projectile  and  sabot  flow  fields. 

An  extensive  experimental  program  [4,5]  was  carried  out  in  1981  to  study  the  aerodynamics  of 
sabot  discard.  Surface  pressures  were  obtained  experimentally  on  a  generic  configuration  consisting 
of  a  projectile  and  three  sabot  components.  The  experimental  test  results  showed  regions  of 
shock/boundary-layer  interactions  and  separated  flow  regions.  Recently,  anodier  experimental  study 
[6]  was  conducted  that  included  detailed  pressure  measurements  on  a  generic  cone-cylinder 
projectile  as  well  as  a  full-scale  model  of  die  C-76  projectile.  The  generic  configuration  was  chosen 
to  provide  calibration  for  computational  fluid  dynamics  (CFD)  modeling  efforts.  As  pointed  out 
earlier,  the  projectile  and  sabot  interacting  flow  field  is  quite  complex  with  3-D  shock/boundary- 
layer  interactions  and  regions  of  separated  flow.  The  CFD  modeling  technique  is  an  emerging  tool 
which  can  account  for  these  3-D  interactions.  Recent  p^rs  [6-8]  show  a  rise  in  the  use  of  CFD  to 
accurately  predict  such  flow  fields.  Often,  the  computed  results  are  found  to  be  in  good  agreement 
with  the  experimental  data.  These  computational  sfeidies  have  provided  enhanced  understanding  of 
the  complex  interacting  flow  fields;  however,  they  used  a  steady-state  approach.  Research  presented 
in  this  report  emphasizes  the  use  of  an  advanced  CFD  capability  that  can  compute  both  steady-state 
and  time-dependent  sabot  discard  aerodynamics. 

The  advanced  CFD  capability  used  here  solves  the  Navier-Stokes  equations  [9]  and  mcorporates 
the  Chimera  technique  [10-13].  The  Chimera  technique  involves  gloating  independent  grids  about 
each  component  and  then  oversetting  them  onto  a  base  grid  to  form  the  complete  model.  With  this 
composite  overset  grid  approach,  it  is  possible  to  use  different  grid  topologies  for  the  projectile  and 
the  sabot  components,  respectively.  A  complete  model  of  the  multibody  system  is  thus  made,  and 
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the  sabot  discard  aerodynamics  can  then  be  determined.  Numerical  flow  field  computations  have 
been  made  for  the  projectile  and  sabot  multibody  system  at  a  supersonic  speed  for  symmetric  sabot 
discard.  Computed  results  have  been  compared  with  the  experimental  data  [6]  obtained  for  the  same 
configuration  and  flow  conditions. 

2.  Solution  Technique 


2.1  Governing  Equations.  The  complete  set  of  S-D,  time-dependent,  generalized  geometiy, 
Reynolds-averaged,  thin-layer  Navier-Stokes  equations  is  solved  numerically  to  obtain  a  solution 
to  this  problem  and  can  be  written  in  general  spatial  coordinates  q,  and  C  as  follows  [14]: 

+  a^F  +  a^6  +  a^ft  =  Re-’a^§,  d) 


hi  Equation  1,  q  contains  the  dependent  variables:  density,  three  velocity  components,  and  energy. 
The  thin-layer  approximation  is  used  here,  and  the  viscous  terms  involving  velocity  gradients  in 
both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous  terms  are  retained 
in  the  normal  direction,  C>  and  are  collected  into  the  vector  S.  These  viscous  terms  are  used 
everywhere. 

2.2  Numnical  Technique.  The  implicit,  approximately  factored  scheme  for  the  thin-layer 
Navi^-Stokes  equations  using  central  differencing  in  the  q  and  C  directions  and  upwinding  in  ^  is 
written  in  the  following  form  [14]: 

[l  .  h6|(A')‘  .  hajC”  -  hRe-'8jJ-'lCi"J  -  D,  |j] 

X  [I +  h6/(A‘)“  +  ha,fl‘  -  D,|,)A0“ 

=  -at(aj[(f  •)"  -  f.*] .  a/[(F')»  -  t.-j  +  a,(G"  -  d.) 

+  aj(ft"  -  ft.)  -  Re-'5j(S"  -  §.))  -  D,(0"  -  0.).  (2) 
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where  /i  =  Af  or  (AO/2  and  the  ftee-stream  base  solution  is  used.  Here,  6  is  typically  a  three-point 
second-order  accurate  central  difference  operator,  6  is  a  midpoint  operator  used  with  the  viscous 
terms,  and  the  operators  and  6^  are  backward  and  forward  three-point  difference  op^tors.  The 
flux  P  has  been  eigensplit,  and  the  matrices  A,  A,  and  lA  result  from  local  linearization  of  the 
fluxes  about  the  previous  time  level.  Here,  J  denotes  the  Jacobian  of  the  coordinate  transformation. 
Dissipation  operators  and  D,-  are  used  in  the  central  space  differencing  directions. 

2.3  Chimera  Composite  Grid  Technique.  The  Chimera  overset  grid  scheme  is  a  domain 
decomposition  approach  where  a  full  configuration  is  meshed  using  a  collection  of  independent 
overset  grids.  This  allows  each  component  of  the  configuration  to  be  gridded  separately  and  overset 
into  a  main  grid.  Overset  grids  are  not  required  to  join  in  any  special  way.  Usually,  a  major  grid 
covers  the  entire  domain  or  a  grid  is  generated  about  a  dominant  body  section.  Minor  grids  are 
generated  about  the  rest  of  the  bodies  or  sections.  Because  each  component  grid  is  generated 
independently,  portions  of  one  grid  may  be  found  to  lie  within  a  solid  boundary  contained  within 
another  grid.  Such  points  lie  outside  the  computational  domain  and  are  excluded  from  the  solution 
process.  Equation  2  has  been  modified  for  Chimera  overset  grids  by  the  introduction  of  the  flag  i^ 
to  achieve  just  that.  This  i^  array  accommodates  the  possibility  of  having  arbitrary  holes  in  the  grid. 
The  i|,  array  is  defined  such  that  i,,  =>  1  at  normal  grid  points  and  i -  0  at  hole  points.  Thus,  when 
ib  “  1,  Equation  2  becomes  the  standard  scheme.  But,  when  i^  -  0,  the  algorithm  reduces  to  A(^"  •= 
0  or  ,  leaving  ^  unchanged  at  hole  points.  The  set  of  grid  points  that  form  the  border 

between  the  hole  points  and  the  normal  field  points  are  called  intergrid  boundary  points.  These 
points  are  updated  by  interpolating  the  solution  from  the  overset  grid  that  created  the  hole.  Values 
of  the  ib  array  and  the  interpolation  coefficients  needed  for  this  update  are  provided  by  a  separate 
algorithm  [10].  The  Chimera  procedure  reduces  a  complex  problem  into  a  number  of  simpler 
subproblems.  Computations  are  performed  on  each  grid  separately.  The  grids  are  developed  to  use 
the  available  core  memory  one  grid  at  a  time.  The  remaining  grids  are  stored  on  an  external  disk 
storage  device  such  as  the  solid-state  disk  (SSD)  device  of  the  Cray  Y-MP  computer.  A  major  part 
of  the  Chimera  overset  grid  approach  is  the  information  transfer  from  one  grid  into  another  by  means 
of  the  intergrid  boundary  points. 
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2  A  Boundary  Conditions.  For  simplicity,  most  of  the  boundary  conditions  have  been  imposed 
explicitly  [9].  An  adiabatic  wall  boundary  condition  is  used  on  the  body  surface,  and  the  no-slip 
boundary  condition  is  used  at  die  wall.  The  pressure  at  the  wall  is  calculated  by  solving  a  combined 
momentum  equation.  Free-stream  boundary  conditions  are  used  at  the  inflow  boundary  as  well  as 
at  the  outer  boundary.  A  symmetry  boundary  condition  is  imposed  at  the  circumferential  edges  of 
the  grid,  while  a  simple  extrapolation  is  used  at  the  downstream  boundary.  A  combination  of 
symmetry  and  extrapolation  boundary  condition  is  used  at  the  carter  line  (axis).  Since  the  free- 
stream  flow  is  supersonic,  a  nonreflection  boundary  condition  is  used  at  the  outer  boundary.  The 
outer  boundary  of  the  sabot  grid  completely  lies  within  the  background  projectile  grid  and,  thus,  gets 
its  flow  field  information  interpolated  from  the  projectile  grid. 

3.  Model  Geometry  and  Computational  Grid 


An  advantage  of  the  Chimera  technique  is  that  it  allows  computational  grids  to  be  obtained  for 
each  body  component  separately  and,  thus,  makes  the  grid  generation  process  easier.  Figure  2  shows 
a  computational  grid  for  the  complete  model,  including  the  projectile  and  sabot.  This  grid 
corresponds  to  the  sabot  angle  of  attack  of  25°.  The  projectile  grid  consists  of  two  zones  (zone  1 
and  zone  2)  that  include  a  small  zone  (zone  1)  in  front  of  the  projectile  nose.  Each  of  these  two  zones 
is  a  rectangular  grid.  The  grid  around  the  sabot  also  consists  of  two  zones  (zone  3  and  zone  4)  and 
was  obtained  using  an  0-topology  and  a  rectangular  topology,  respectively.  Figure  2  also  shows  the 
sabot  grid  (zone  3).  The  sabot  grids  were  individually  generated  and  then  overset  as  shown  in  this 
figure  to  form  the  complete  grid  system.  The  computational  grids  shown  here  correspond  to  the  pitch 
plane.  The  zone  2  projectile  grid  serves  as  the  main  background  grid  for  the  computation.  The  grid 
in  zone  3  is  an  O-grid  around  the  sabot  petal.  The  grid  in  zone  4  sits  along  the  edge  of  the  sabot 
petal  (not  shown  in  this  figure).  Figure  3  shows  a  computational  grid  for  computations  with  the 
sabot  petal  at  angle  of  attack  of  10° .  Again,  the  same  sabot  grids  are  used  hoe  for  this  run,  and  there 
was  no  need  to  regenerate  new  sabot  grids.  The  same  sabot  grids  were  also  used  for  two  other  cases, 
which  correspond  to  sabot  angles  of  attack  of  5°  and  15°.  In  each  case,  the  dimensions  of  different 
zones  are  as  follows:  zone  1,  16x32x80;  zone  2,  86x32x80;  zone  3,  72x32x30;  and  zone  4, 
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60x20x40.  The  entire  grid  system  consisted  of  378,240  grid  points.  Note  that  the  grid  setup  allows 
computation  of  the  base  region  flow  field  of  the  sabots.  Grid  points  are  clustered  near  the  projectile 
and  the  sabot  surfaces  to  capture  the  viscous  boundary  layers.  No  attempt  has  been  made  to  adapt 
the  computational  grids  to  gradients  in  the  flow  field  vaiiables. 


Figure  2.  Computational  Grid,  Sabot  Angle  of  Attack  =  25°. 


Figure  3.  Computational  Grid,  Sabot  Angle  of  Attack  =  10°. 

4.  Results 

Steady-state  numerical  calculations  have  been  performed  to  numerically  simulate  the  projectile 
and  the  sabot  system.  All  computations  have  been  run  at  =  4.0  and  for  the  same  test  conditions 
corresponding  to  the  Canadian  experiments  [6].  Computational  modeling  is  restricted  to  the 
symmetric  sabot  discard.  Here,  the  projectile  is  at  zero  degrees  angle  of  attack  and  three  sabots  are 
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discarded  symmetrically  following  the  same  radial  trajectory  away  from  the  projectile.  Figure  4 
shows  a  schematic  diagram  of  the  projectile  and  the  sabot  system.  Since  symmetric  sabot  discard  is 
of  interest  here,  the  computational  domain  consists  of  a  60®  segment,  as  shown  in  Figure  4.  Also 
shown  is  the  sabot  grid,  which  is  entirely  contained  in  the  background  projectile  grid.  Because  of 
symmetry,  the  requirements  for  grid  sizes,  computer  resources  such  as  computer  memory,  and  run 
time  are  reduced. 


Figure  4.  Schematic  Diagram  Showing  the  Computationai  Domain 


Initially,  a  converged  result  for  the  sabot  at  25°  angle  of  attack  was  obtained.  Figure  5  shows 
the  computed  Mach  contours  for  this  case.  It  shows  the  computed  Mach  contours  for  a 
circumferential  cross-sectional  plane  which  cuts  through  the  projectile  and  the  sabots.  Although 
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Figure  5.  Shock  Interactions  in  the  Cross-Sectionai  Piane  (Projectiie  and  Sabots). 


Figure  6.  Computed  Pressure  Contours,  M  =  4.0,  Sabot  Angle  of  Attack  =  25®,  Laminar  and 
Turbulent  Solutions. 
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computations  are  pafonned  only  for  a  60“  sector  due  to  symmetry,  the  computed  results  are  shown 
here  for  the  entire  projectile  and  the  sabot  system.  This  figure  clearly  shows  the  shock  interactions 
associated  with  the  projectile  and  the  sabot  s)«tem.  Figure  6  shows  the  computed  pressure  contours 
for  the  25“  sabot  angle  of  attack  case  in  the  pitch  plane.  Both  laminar  and  turbulent  computations 
were  performed  and  only  small  differences  were  observed  between  these  results,  especially  near  the 
sabot  shock  impingement  point  on  the  projectile.  The  turbulent  calculations  show  slightly  larger 
regions  of  shock-boundary  layer  interactions.  The  wind  tunnel  data  for  this  case,  however,  was  later 
found  to  have  suifraed  from  the  Wall  blockage  effects,  and  thus,  numerical  results  for  this  case  were 
not  compared  with  the  data.  It  shows  an  oblique  shock  wave  emanating  from  the  nose  of  the 
projectile  and  a  detached  bow  shock  in  front  of  the  sabot  petal.  The  interactions  of  the  projectile  and 
the  sabot  flow  fields  can  be  seen  clearly.  Downstream  of  the  shock-interaction  point,  a  region  of 
high  pressure  and  low  velocity  exists.  Based  on  the  shock  strengths,  this  interacting  flow  region  may 
include  regions  of  flow  separation  in  addition  to  shock-boundary  layer  interactions.  The  computed 
result  shown  corresponds  to  the  symmetry  plane.  A  nonreflection  outer  boundary  is  used,  which 
allows  the  outer  boundary  of  the  background  computational  grids  to  be  placed  near  the  sabots. 

The  next  set  of  computed  results  corresponds  to  the  B5,  BIO,  and  B1 5  cases.  These  three  cases 
correspond  to  sabot  angles  of  attack  of  5“,  10“,  and  1S“.  The  projectile  is  at  zero  degrees  angle  of 
attack  for  all  these  cases.  As  stated  earlier,  the  background  grid  for  the  projectile  remains  the  same. 
The  sabot  grids  again  are  the  same  but  have  been  moved  to  the  new  positions  and  orientations. 
Figure  7  shows  the  pressure  contours  for  the  projectile  and  sabot  in  the  symmetry  plane.  Computed 
results  here  have  been  obtained  for  the  turbulent  flow  condition  using  an  algebraic  turbulence  model. 
Similar  to  the  25“  sabot  angle  of  attack  case,  a  nonreflection  outer  boundary  is  used,  which  allows 
the  outer  boundary  of  the  background  computational  grids  to  be  placed  near  the  sabots.  This  figure 
clearly  shows  the  interactions  of  the  projectile  and  the  sabot  flow  fields  occurring  at  different 
longitudinal  locations  along  the  projectile.  The  computed  pressure  contours  show  the  sabot  shock 
impinging  and  reflecting  off  the  projectile  surface.  This  shock  impingement  results  in  a  higher 
pressure  region  on  die  projectile  surface  just  downstream  of  the  impingement  point.  A  high  pressure 
region  can  also  be  observed  behind  the  sabot  shock.  As  expected,  the  flow  behind  the  base  region 
of  the  sabot  is  a  low-pressure  region.  As  the  sabot  angle  of  attack  is  increased,  the  sabot  shock 
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impingement  point  on  the  projectile  is  moved  further  downstream.  For  the  5°  sabot  angle  of  attack, 
the  sabot  shock  impinges  on  the  projectile,  reflects  from  the  projectile  surface,  and  impinges  back  on 
the  sabot.  The  reflected  shock  from  the  projectile  surface  is  seen  to  just  miss  the  base  of  the  sabot 
for  the  10°  sabot  angle  of  attack  case  and  is  even  further  away  from  the  sabot  base  at  15°  sabot  angle 
of  attack.  The  flow  field  in  the  base  region  of  the  sabot  is  also  seen  to  change  considerably  with  an 
increase  in  sabot  angle  of  attack. 


Figure  7.  Computed  Pressure  Contours  at  Different  Sabot  Angles  of  Attack,  5°,  10°,  and  15°, 
M  =  4.0. 
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Figures  8  and  9  show  the  computed  surface  pressures,  p/p„  for  the  sabot  and  the  projectile, 
respectively.  Here,  X/D  ■  0.0  corresponds  to  the  nose  of  the  projectile.  These  computed  surface 
pressures  correspond  to  the  pitch  plane  and  are  compared  with  the  experimental  data  [6].  The 
computed  pressures  on  the  bottom  surface  of  the  sabot  are  shown  in  Figure  8  and  are  generally  in 
good  agreement  with  the  experimental  data.  Some  discrepancies  do  exist  in  the  comparison  of  sabot 
surface  pressure  for  the  5”  angle  of  attack  case.  Due  to  close  proximity  of  ttie  sabot  to  the  projectile, 
the  flow  field  is,  as  expected,  more  complicated  and  includes  complex  shock>shock  and  shock- 
boundary  layer  interactions.  Accurate  computation  of  die  resulting  flow  field  is  thus  more  difficult. 
Grid  clustering/alignment  in  the  nose  region  of  the  sabot  may  be  needed  to  improve  the  accuracy  of 
the  numerical  solution  for  this  case.  Similar  experience  has  been  noted  by  other  researchers  [6].  For 
X/D  >  7,  the  agreement  of  the  computed  surface  pressures  with  the  data  is  very  good  even  for  the 
5“  angle  of  attack  case.  For  this  case,  there  is  a  secondary  pressure  rise  near  X/D  =  11 ,  which  results 
from  the  reflected  shock  impinging  on  the  underside  of  the  sabot.  This  predicted  pressure  rise 
matches  very  well  with  the  experimental  data.  As  the  sabot  angle  of  attack  is  increased,  the 
secondary  pressure  rise  is  reduced  until  it  is  eliminated  at  the  15®  sabot  angle  of  attack.  The 
agreement  of  the  computed  sabot  surface  pressures  with  the  data  is  good  at  angle  of  attack  of  10® 
and  15®.  Figure  9  shows  the  surface  pressure  distributions  on  the  projectile  in  the  pitch  plane. 
Computed  results  are  shown  in  solid  line  and  are  compared  with  the  experimental  data  shown  in  dark 
circles  for  5®,  10®,  and  15®  degree  sabot  angles  of  attack.  As  seen  in  this  figure,  the  surface  pressure 
is  almost  constant  on  the  nose,  which  is  followed  by  a  pressure  drop  at  the  cone-cylind^  junction. 
This  computed  pressure  drop  at  the  cone-cylinder  junction  agrees  well  with  the  data  at  the  15®  sabot 
angle  of  attack;  however,  at  lower  angles  of  attack,  die  agreement  is  not  so  good.  The  predicted  flow 
on  the  nose  of  die  projectile  corresponds  to  an  undisturbed  flow  upstream  of  the  shock  impingement 
point.  Clearly,  the  numerical  results  do  not  show  the  same  extent  of  shock-boundary  layer 
interactions  observed  experimentally.  Similar  results  are  also  seen  in  the  DREV  CFD 
predictions  [6].  A  large  pressure  increase  due  to  shock  wave  impinging  on  the  projectile  surface  is 
seen  in  both  computed  and  experimental  data.  The  locations  of  the  pressine  peaks  have  been 
predicted  correcdy  and  agree  well  with  the  data.  The  magnitudes  of  the  peak  are,  however,  slightly 
underpredicted.  Additionally,  a  secondary  small  peak  is  observed  in  the  experimentally  obtained 
surface  pressure  near  X/D  =  1 1  for  the  5®  case,  which  is  not  se«i  in  the  computed  surface  pressures. 
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It  may  be  possible  to  improve  the  accuracy  of  the  computed  results  eith^  through  grid  refinement 
and/or  use  of  advanced  turbulence  modeling,  but  further  analysis  has  been  limited  due  to  time 
constraints.  ;  * 


SABOT  B5  PROiECTILCBS 

K 


SABOT  BIO  MOJKTILEOtO 


SABOT  BtS  PROJECTILE  B» 


Figure  8.  Sabot  Surface  Pressure 


Figure  9.  Projectile  Surface  Pressure 


Distributions. 


Distributions. 
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5.  Concluding  Remarks 


A  computational  study  was  undertaken  to  compute  the  3-D  flow  fields  for  a  multibody  system 
consisting  of  a  projectile  and  sabots.  Flow  computations  were  performed  at  a  supersonic  Mach 
number,  M„  -  4.0  and  a  -  0.0®  using  a  3-D  unsteady  Navier-Stokes  code  and  Chimera  composite 
grid  discretization  technique.  Overset  body  conforming  grids  were  used  to  individually  model  the 
projectile  and  the  sabot  components.  Computed  r^ults  have  been  obtained  for  sabot  angles  of  attack 
of  5®,  10®,  15®  and  25®.  Computed  results  show  the  qualitative  features  of  the  complex  shock 
interaction  flow  field  for  the  projectile  and  the  sabots.  Both  laminar  and  turbulent  computed  results 
have  been  obtained  for  the  25®  case,  and  the  computed  results  do  not  show  appreciable  change  in 
the  surface  pressures.  Computed  results  for  this  case  have  not  been  compared  widi  the  experimental 
results  due  to  blockage  effects  encountered  in  the  tests.  Computed  results  for  the  other  three  sabot 
positions  are  compared  with  the  experimental  data  obtained  at  DREV,  Canada,  for  the  same 
configuration  and  conditions  and  are  generally  found  to  be  in  good  agreement  with  the  data.  In  some 
cases,  discrepancies  exist  between  the  computed  surface  pressures  and  experimental  data.  Grid 
refinement  and  use  of  advanced  turbulence  modeling  may  be  heeded  to  furdier  improve  the  accuracy 
of  the  computed  results.  FUture  study  will  include  modeling  of  asymmetric  sabot  discard,  which  will 
require  full  3-D  computations  and  large  computing  resources. 


Intentionally  l^  blank. 


6.  References 


1 .  Schmidt,  E.  M.,  and  D.  D.  Shear.  "Aerodynamic  Interference  During  Sabot  Discard."  AIAA 
Journal  of  Spacecraft  and  Rockets,  vol.  15,  no.  3,  May-June  1978,  pp.  162-167. 

2.  Crimi,  R,  and  D.  Siegelman.  “Analysis  of  Mechanical  and  Gasdynamic  Loadings  During  Sabot 
Discard  From  Gun-Launched  Projectiles."  ARBRL-CR-341,  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD,  June  1977. 

3.  Siegelman,  D.,  J.  Wang,  and  P.  Crimi.  “Computation  of  Sabot  Discard."  ARBRl^CR-505,  U.S. 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  February  1983. 

4.  Schmidt,  E.  M.  “Wind-Tunnel  Measurements  of  Sabot  Discard  Aerodynamics."  AIAA  Journal 
of  Spacecraft  and  Rockets,  vol.  18,  no.  3,  May-June  1981,  pp.  235-242. 

5.  Schmidt,  E.  M.,  and  P.  Plostins.  “Aerodynamics  of  Asymmetric  Sabot  Discard."  ARBRI^ 
MR-03281,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD, 

June  1983. 

6.  Lesage,  F.,  and  B.  Girard.  “Wind  Tiinnel  and  CFD  hivestigation  of  Aerodynamic  hiteractions 
During  Sabot  Separation."  AIAA  Paper  No.  96-0193,  January  1996. 

7.  Nusca,  M.  “Computational  Fluid  Dynamics  Application  to  the  Aerodynamics  of  Symmetric 
Sabot  Discard.”  AIAA  Paper  No.  90-3096,  August  1990. 

8.  Nusca,  M-  “Numerical  Simulation  of  Sabot  Discard  Aerodynamics.”  AIAA  Paper  No.  91-3255, 
September  1991. 

9.  Sahu,  J.  "Numerical  Computations  of  Transonic  Critical  Aerodynamic  Behavior."  AIAA 
Journal,  vol.  28,  no.  5,  pp.  807-816,  May  1990. 

10.  Steg«:,  J.  L.,  F.  C.  Dougherty,  and  J.  A.  Benek.  "A  Chimera  Grid  Scheme."  Advances  in  Grid 
Generation,  edited  by  K.  N.  Ghia  and  U.  Ghia,  ASME  FED-5,  June  1983. 

11.  Benek,  J.  A.,  T.  L.  Donegan,  and  N.  E.  Suhs.  'Extended  Chimera  Grid  Embedding  Scheme 
With  Application  to  Viscous  Flows."  AIAA  Paper  No.  87-1 126-CP,  1987. 

1 2.  Meakin,  R.  L.  "Computations  of  the  Unsteady  Flow  About  a  Generic  Wing/Pylon/Finned-Store 
Configuration."  AIAA  92-4568-CP,  August  1992. 

13.  Sahu,  J.,  and  C.  J.  Nietubicz.  "Application  of  Chimera  Technique  to  Projectiles  in  Relative 
Motion."  ARL-TR-590,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD, 
October  1994  (also  see  AMA  Journal  of  Spacecraft  and  Rockets,  vol.  32,  no.  5,  Sep-Oct  1995). 


15 


14.  Steger,  J.  L.,  S.  X.  Ying,  and  L.  B.  Schiff.  "A  Partially  Flux-Split  Algorithm  for  Numerical 
Simulation  of  Compressible  inviscid  and  Viscous  Flows."  Proceedings  of  the  Woritshop  on 
CFD,  Institute  of  Nonlinear  Sciences,  Universi^  of  California,  Davis,  CA,  1986.. 


NO.  OF 

COPES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  DEFENSE  TECHNICAL 
INFORMATION  CENTER 
DTIC  DDA 

8725  JOHN  J  KINGMAN  RD 
STE0944 

FT  BELVOIR  VA  22060-6218 

1  HQDA 

DAMOFDQ 
DENNIS  SCHMIDT 
400  ARMY  PENTAGON 
WASHINGTON  DC  20310-0460 

1  CECOM 

SP  &  TRRSTRL  COMMCTN  DIV 

AMSELRDSTMCM 

HSOICHER 

FT  MONMOUTH  NJ  07703-5203 

1  PRIN  DPTY  FOR  TCHNIGY  HQ 

USARMYMATCOM 
AMCDCGT 
M  FISETTE 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  PRIN  DPTY  FOR  ACQUSTN  HQS 
USARMYMATCOM 
AMCDCGA 
D  ADAMS 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  DPTY  CG  FOR  RDE  HQS 

USARMYMATCOM 
AMCRD 

BG  BEAUCHAMP 

5001  EISENHOWER  AVE 

ALEXANDRIA  VA  22333-0001 

1  DPTY  ASSIST  SCY  FOR  R&T 

SARDTT  TKILUON 
THE  PENTAGON 
WASHINGTON  DC  20310-0103 

1  OSD 

OUSDCA&iyODDDR&ECR) 

JLUPO 

THE  PENTAGON 
WASHINGTON  DC  20301-7100 


1  INST  FOR  ADVNCD  TCHNIGY 

THE  UNIV  OF  TEXAS  AT  AUSTIN 
PO  BOX  202797 
AUSTIN  TX  78720-2797 

1  USAASA 

MOASAI  WPARRON 
9325  GUNSTON  RD  STE  N319 
FT  BELVOIR  VA  22060-5582 

1  CECOM 

PM  GPS  COLS  YOUNG 
FT  MONMOUTH  NJ  07703 

1  GPS  JOINT  PROG  OFC  DIR 
COL  J  CLAY 

2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA  90245-5500 

1  ELECTRONIC  SYS  DIV  DIR 
CECOM  RDEC 
JNIEMELA 

FT  MONMOUTH  NJ  07703 

3  DARPA 
L  STOTTS 
JPENNELLA 
BKASPAR 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

1  SPCL  ASST  TO  WING  CMNDR 

50SW/CCX 

CAPT  PH  BERNSTEIN 
300  O'MALLEY  AVE  STE  20 
FALCON  AFB  CO  80912-3020 

1  USAFSMaCED 
DMA/JPO 
MISON 

2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA 
90245-5500 

1  US  MIIiTARY  ACADEMY 

MATH  Sa  CTR  OF  EXCELLENCE 
DEPT  OF  MATHEMATICAL  SCI 
MDN  A  MAJ  DON  ENGEN 
THAYER  HALL 
WEST  POINT  NY  10996-1786 


17 


NO.  OF 

COPIES  ORGANIZATION 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLCSALTP 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1145 

1  DIRECrOR 

US  ARMY  RESEARCH  LAB 
AMSRLCSALTA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRLCILL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 


ABERDEEN  PROVING  GROUND 

4  DIRUSARL 

AMSRLCILP(305) 


NO.  OF 

COPIES  ORGANIZATION 


2  USAF  WRIGHT  AERONAUTICAL 

LABORATORIES 
ATTN  AFWALFIMG 
DRJSHANG 
MRNESCAGGS 
WPAFB  OH  45433-6553 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CNTR 
ATTN  CODE  B40  DR  W  YANTA  ' 
DAHLGREN,  VA  22448-5100 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CNTR 
ATTN  CODE420DR  AWARDLAW 
INDIAN  HEAD,  MD  20640-5035 

4  DIRECTOR 

NATIONAL  AERONAUTICS  AND 
SPACE  ADMINISTRATION 
LANGLEY  RESEARCH  CENTER 
ATTN  TECH  LIBRARY 
MRDMBUSHNELL 
DRMJHEMSCH 
DRJSOUTH 
LANGLEY  STATION 
HAMPTON  VA  23665 

2  ARPA 

ATTN  DRPKEMMEY 
DR  JAMES  RICHARDSON 
3701  NORTH  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

7  DIRECTOR 
NASA 

AMES  RESEARCH  CENTER 
ATTN  MS  227  8  L  SCHIFF 
MS  258  IT  HOLST 
MS2581DCHAUSSEE 
MS2581MRAI 
MS2581PKUTLER 
MS258  1PBUNING 
MS  258  1  B  MEAKIN 
MOFFETT  FIELD  CA  94035 


NO.  OF 

COPIES  ORGANIZATIQN 


2  USMA 

DEPT  OF  MECHANICS 
ATTN  LTC  ANDREW  L  DULL 
M  COSTELLO 
WEST  POINT  NY  10996 

7  COMMANDER 
USARMYARDEC 
ATTN  SMCARAETA 
RDEKLEINE 
CNG 

R  BOTTICELLI 
H  HUDGINS 
JGRAU 
SKAHN 
W  KOENIG 

PICATTNNY  ARSENAL  NJ  07806-5001 

1  COMMANDER 
USARMYARDEC 
ATTN  SMCARCCHV 
PAUL  VALENTI 

PICATINNY  ARSENAL  NJ  07806-5001 

1  COMMANDER 
USARMYARDEC 
ATTN  SFAEFASSD 
MIKE  DEVINE 

PICATINNY  ARSENAL  NJ  07806-5001 

1  COMMANDER 

US  NAVAL  SURFACE  WEAPONS  CTR 
ATTN  DR  F  MOORE 
DAHLGREN  VA  22448 

2  UNIV  OF  CALIFORNIA  DAVIS 
DEPT  OF  MECHANICAL  ENGG 
ATTN  PROF  HA  DWYER 
PROF  M  HAFEZ 

DAVIS  CA  95616 

1  AEROJET  ELECTRONICS  PLANT 
ATTN  DANIEL  W  PILLASCH 
B170  DEPT  5311 
PO  BOX  296 

1 100  WEST  HOLLYVALE  STREET 
AZUSA  CA  91702 


19 


NO.  OF 

corns  QRGANIZATIQN 

3  SCIENCE  AND  TECHNOLOGY  INC 
4001  NORTH  FAIRFAX  DR  NO  700 
ATTN  DRALANGLASSER 
MR  BRUCE  LOHMAN 
MRDAVEMAURIZI 
ARLINGTON  VA  22203-1618 

3  AIR  FORCE  ARMAMENT  LAB 
ATTN  AFATiyFXA 
STEPHEN  C  KORN 
BRUCE  SIMPSON 
DAVEBELK 

EGLINAFB  FL  32542-5434 

1  MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
ATTN  TECH  LIBRARY 
77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 

1  GRUMANN  AEROSPACE  CORPORATION 
AEROPHYSICS  RESEARCH  DEPT 
ATTN  DRREMEU4IK 
BETHPAGE  NY  11714 

2  MICRO  CRAFT  INC 
ATTN  DR  JOHN  BENEK 
NORMAN  SUHS 

207  BIG  SPRINGS  AVE 
TULLAHOMA  TN  37388-0370 

1  LOS  ALAMOS  NATIONAL  LAB 

ATTN  MR  BILL  HOGAN 
MS  G770 

LOS  ALAMOS  NM  87545 

3  DIRECTOR 

SANDIA  NATIONAL  LABORATORIES 
’  ATTN  DIV1554DRWOBERKAMPF 
DIV  1554  DR  F  BLOTTNER 
DIV  1636  DRW  WOLFE 
ALBUQUERQUE  NM  87185 

1  NAVAL  AIR  WARFARE  CENTER 

ATTN  DAVID  FINDLAY 
MS  3  BLDG  2187 
PATUXENT  RIVER  MD  20670 


NO.  OF 

COPIES  ORGANIZATION  • 

1  METACOMP  TECHNOLOGIES  INC 
ATTN  S  R  CHAKRA  VARTHY 

650  HAMPSHIRE  ROAD 
SUITE  200 

WESTLAKE  VILLAGE  CA  91361-2510 

2  ROCKWELL  SCIENCE  CENTER 
ATTN  S  V  RAMAKRISHNAN 

V  V  SHANKAR 
1049  CAMINO  DOS  RIOS 
'  THOUSAND  OAKS  CA  91360 

1  ADVANCED  TECHNOLOGY  CTR 

ARVIN/CALSPAN  , 
AERODYNAMICS  RESEARCH  DEPT 
ATTN  DR  MS  HOLDEN 
PO  BOX  400 
BUFFALO  NY  14225 

1  PENNSYLVANIA  STATE  UNIV 
DEPT  OF  AEROSPACE  ENGG 
ATTN  DRGSDULKRAVICH 
UNIVERSITY  PARK  PA  16802 

1  UNIV  OF  ILLINOIS  AT  URBANA 

CHAMPAIGN 

DEPT  OF  MECHANICAL  AND 
INDUSTRIAL  ENGINEERING 
ATTN  DR  JC  DUTTON 
URBANA  IL  61801 

1  UNIVERSITY  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGG 
ATTN  DR  J  D  ANDERSON  JR 
COLLEGE  PARK  MD  20742 

1  UNIVERSITY  OF  NOTRE  DAME 

DEPT  OF  AERONAUTICAL  AND 
MECHANICAL  ENGINEERING 
ATTN  PROF  TJ  MUELLER 
NOTRE  DAME  IN  46556 

1  UNIVERSITY  OF  TEXAS 

DEPT  OF  AEROSPACE  ENGG 
MECHANICS 
ATTN  DRDSDOLLING 
AUSTIN  TX  78712-1055 


20 


1 


UNIVERSITY  OF  DELAWARE 
DEPT  OF  MECHANICAL  ENGG 
ATTN  DRJOHNMEAKIN 
NEWARK  DE  19716 

1  UNIVERSITY  OF  FLORIDA 
DEPT  OF  ENGG  SCIENCES 
COLLEGE  OF  ENGINEERING 
ATTN  PROF  CC  HSU 
GAINESVEJLE  FL  326i  1 


ABERDEEN  PROVING  GROUND 

27  DIR  USARL 

ATTN  AMSRL-WM-P 
A HORST 
E  SCHMIDT 
AMSRL-WM-PB 
PPLOSTINS 
DLYON 
MBUNDY 
KFANSLER 
EFERRY 
B  GUIDOS 
KHEAVEY 
HEDGE 
V  OSKAY 
AMIKHAIL 
JSAHU 
PWEINACHT 
AMSRL-ST,  J  ROCCHIO 
AMSRL-WM-PD,  B  BURNS 
AMSRL-WM-PA 
GKELLER 
MNUSCA 

AMSRL-WM-PC,  B  FORCH 
AMSRL-WM-W,  C  MURPHY 
AMSRL-WM-WB,  W  D’AMICO 
AMSRL-WM-TB,  R  LOTTERO 
AMSRL-CI-H 
CNIETUBICZ 
AMSRL-CI-HC 
P  COLLINS 
DHISLEY 
DPRESSEL 
WSTUREK 

2  CDR  ARDEC 

ATTN  FIRINGTABLES,  BLDG  120 
RLIESKE 
REITMILLER 


21 


Intentionally  left  blank. 


22 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  Na  0704^188 


I  lTi  7 1 ;  r^Ti  rnr; 


1  itour  per  raaponM.  inciitoing  m*  tima  ler  rMmflng  mnrvevonn,  wemnm^  ensung  onv  wurpw. 
gotlwlng  and  maliitalning  Vm  d«a  naadod.  and  eomplaang  and  revlawtog  the  eollaetton  af  formation.  Sand  commanfs  raganltng  tNs  fauiden  aatimata  or  any  ottiar  aapaet  of  this 
collaetlon  of  Infonnailon,  inefudng  auggeattona  lof  radudRg  iMs  buntan,  to  WaatiEngton  Headquarters  Sarvtees,  OtrectoraCa  for  tmtonnattM  Oparatiorta  and  Reporiai  121S  JaMafaon 
Oavte  Hkihwaw.  Suita  1204.  AfBnotcn.  VA  22202-4302.  end  to  the  Offtce  of  Management  and  Budget  Pepenwork  Beductlen  Protect(OT04»OlB».  Waablnoton.  DC  aosoa. _ 


1.  AGENCY  USE  ONLY  (Lews  bfanl^  2.  REPORT  DATE 

October  1997 


3.  REPORT  TYPE  AND  DATES  COVERED 

Pinal,  Oct  95  -  Sep  96 


Navier-Stokea  Computations  of  Sabot  Discard  Using  Chimera  Sdieme 


5.  FUNDING  NUMBERS 

1L161102AH43 


E.  Feny,  J.  Sahu,  and  K.  Heaver 


RGANJZATI 

U.S.  Anny  Research  Laboratoiy 

ATTN:  AMSRI^WM  PB 

Aberdeen  Proving  Ground,  MD  21005-5066 


a.  PERFORMING  ORGANIZATIO 
REPORT  NUMBER 


ARL-TR-1524 


0.  SPONSORINOmONITORINQ  AGENCY  NAMES{S)  AND  ADDR 


10.8PON8ORiN 

AGENCY  REPORT  NUMBER 


11.  SUPPLB4ENTARY 


12a.  DISTRIBUTION/AVAILABIUTY^ATEMENT 
ApjHOved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (MaxftttCfm  200  words) 

Computational  fluid  dynamics  calculations  have  been  performed  for  a  multibody  system  consisting  of  a  main 
projectile  and  three  sabot  components.  Numerical  flow  field  computations  have  been  made  for  various  orientations 
and  locations  of  sabots  using  an  unsteady,  zonal  Navier-Stokes  code  and  the  Chimera  composite  grid  discretizatioa 
technique  at  M.  4.0  and  a  -  0°.  Computational  grids  have  been  obtained  for  foe  projectile  and  sabot 
ind^)endently  and  then  overset  to  form  the  complete  grid  system.  Computed  results  have  been  obtained  for  sabot 
angles  of  attack  of  5, 10,  IS,  and  25°.  Computed  results  show  foe  details  of  foe  expected  flow  field  features 
including  foe  shock  interactions.  Bofo  laminar  and  turbulent  computations  for  the  25°  case  predict  similar  results. 
Computed  results  for  other  sabot  positions  are  compared  with  the  experimental  data  obtained  in  Canada  for  foe 
same  configuration  and  conditions  and  are  generally  found  to  be  in  good  agreement  with  foe  data. 


14.  SUBJECT  TERMS  18.  NUMBER  OF  PAGES 

kinetic  energy  projectiles,  sabot  discard,  Navier-Stokes  equations.  Chimera  composite  25 

grid  technique,  steady-state  numerical  calculations,  shock  interaction,  pressure  contours,  le.  price  code 
surface  pressure  distrifauti<ai _ 


17.  security  CLASSIHCATION  18.  SECURITY  CLASSIFICATION  1».  SECURITY  CLASSIFICATION  20.  UMITATION  OF  ABSTRACT 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED  UL 


7540-01-280-6500  Standard  Form  298  (Rev.  2-89} 

Prescribed  by  ANSI  SM.  239-18  298-102 


23 


Intentionally  left  blank. 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1 .  ARL  Report  Number/ Author  ARL-TR-1 524  fFerrv') _ Date  of  Report  October  1997 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report  will 

be  used.) _ ^ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate.  _ _ _ _ 


6.  General  Commits.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to  organization, 
technical  content,  format,  etc.)  _ _ ^ _ 


Organization 

CURRENT 

ADDRESS 

Name 

E-mail  Name 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 


7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Currrat  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001,APG,MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRLWMPB 

ABERDEEN  PROVING  GROUND  MD  21005’5066 


